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Abstract. We establish a quenched Central Limit Theorem (CLT) for a 
smooth observable of random sequences of iterated linear hyperbolic maps 
on the torus. To this end we also obtain an annealed CLT for the same sys- 
tem. We show that, almost surely, the variance of the quenched system is the 
same as for the annealed system. Our technique is the study of the transfer 
operator on an anisotropic Banach space specifically tailored to use the cone 
condition satisfied by the maps. 



1. INTRODUCTION 

The issue of limit laws in dynamical systems has been widely explored in the last 
decades and it has a clear relevance for physical applications. A prime example of 
a physically relevant system is the study of the statistical behavior of a Lorenz gas 
with randomly distributed obstacles. The case of periodic obstacles is known to 
be ergodic. This follows from the recurrence [19], which in turns follows from the 
CLT, [10, 20], which has been proved in [18] (see also [8] for more refined results 
on these issues). See [15] for more details and for the treatment of some (locally) 
aperiodic cases. On the contrary the random case (albeit one may naively think of 
it as an easier case) stands as a challenge. 

If one considers the simplest possibility (the random position of the obstacles is 
a small i.i.d. perturbation of a periodic configuration) then, by Poincare section, 
one is readily reduced to considering a random sequence of hyperbolic symplcctic 
maps. Yet, such a sequence of maps is not i.i.d. due to the presence of recollisions. 
Recollisions are notoriously a source of serious problems in the study of gases but, 
quite surprisingly, even disregarding the recollision problem (i.e. for the i.i.d. case), 
the problem is poorly understood. 

In this paper we address the easiest setting in which such a situation occurs: 
an i.i.d. sequence of smooth uniformly hyperbolic symplectic maps. To make the 
presentation as clear as possible we will steer away from the full generality in 
which the present results can be obtained (although we will comment on it) and 
we will consider an i.i.d. sequence of linear two dimensional toral automorphisms. 
Exponential decay of correlations has been shown in this setting in [1]. 

For such a model we will show that the time-iV average of any smooth zero 
mean observable has Gaussian fluctuations of order s/N for almost every sequence 
of maps. Moreover, we identify the variance of such Gaussian fluctuations. 
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Similar, but less complete, results are obtained in [2] where the Gaussian na- 
ture of the fluctuations is established for each sequence of maps but neither the 
amplitude nor the variance is investigated. 

The paper is self-contained and organized as following: Section 2 contains the 
precise description of the model we will discuss and states the main results of 
the paper. Such results depend on the understanding of the ergodic properties of 
sequences of maps. These are investigated in Section 3 where the needed ergodic 
properties are related to the spectral properties of transfer operators viewed on 
appropriate Banach spaces in the spirit of the line of research started with [6]. 
Next, in Section 4, we use the above results to establish a CLT averaged over the 
environment for a class of systems larger than the ones at hand but necessary to 
handle the quenched case. The latter is dealt with in Section 5 using an approach 
inspired by works on random walks in random environments; see [11] and references 
therein. 

Convention 1.1. In this paper we will use C to designate a generic constant 
depending only on the choice of the maps {T!;} below. We will use instead C a ^. c ,... 
for constants depending also on the parameters a, 6, c, . . . . Consequently, the actual 
numerical value of such constants may vary from one occurrence to the next. On 
the contrary we will use C\,C2, ■ ■ ■ , to designate constants whose value is fixed 
throughout the paper. 

Acknowledgments. Wc would like to thank Joel Lebowitz for posing the problem 
and Dmitry Dolgopyat for communicating to one of us (CL) reference [2]. MS 
would like to thank the Finnish Cultural Foundation for funding. MS and AA were 
supported in part by NSF DMR-01-279-26 and AFOSR AF 49620-01-1-0154. CL 
would like to thank the Courant Institute where he was visiting when this work 
started. 



2. The model and the results 

Let us consider two 1 matrices {^L}i =0 G SL(2,N) and define the toral automor- 
phisms Tix = AjX mod 1. Let p e [0, 1] and set po = p, p\ = 1 — p. We can then 
introduce the Markov operator Q p : L°°(T 2 ,R) -> i°°(T 2 ,K) defined by 

l 

Qp9( x ) = ^2pi9(Ti(x)). 

i=0 

Such an operator defines a Markov Process. To describe it wc consider the space of 
trajectories J7* := (T 2 ) N endowed with the product topology and, letting (xq, Xi, . . . ) 
be a general element in f2*, we have the obvious dynamics t : f)» — > fi» defined 
by t(.to, xi, . . .) — (xi, . . . ). For each initial measure \i on T 2 , the above Markov 
process defines a Borel probability measure on f2». Let Ep be the expectation 
with respect to such a measure. Then 

Ep^a*)) = J g(t)Kdt) 
E p fl (gfa+i) I x i) = Qpg(xi)- 



In fact, the following would hold almost verbatim also for any larger collection of matrices. 
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The measure P^ is supported on a very small set of trajectories: P^-almost surely 
x i+1 G {T Xi,TiXi}. Thus, if we consider S := {0, l} N Vf°} ; we have 

-P^(U we s{(^o, T Ul xo, T ul2 T u>1 xo, ...)}) = 1. 

In other words we can define the probability space CI = X x T 2 (again equipped 
with the product topology), the map F : CI — > CI defined by 

F(u,x) = (tu,T Ui x), 

where (rw)j = Wj+i, and the measure P P4l — P p x fi, where P p is the Bernoulli 
measure with probability p of having zero. We will denote by E p the expectation 
with respect to P p . Note that if fi is simultaneously To an d Ti invariant, then 
P PiA1 is invariant for the map F. Since the maps are symplectic, this happens for 
the normalized Lebesgue measure m. Let us set P p := P p . m and call Ep p the 
corresponding expectation. Finally, we define the map \I r : Cl — * f2* by 

x) := (x,T Ul x,T U3 T Ul x, . . .). 

It is then easy to verify that T k (^(tu, x)) = ^(F k (oj, x)) for all (u,x) G O, k G N, 
and Ep (/i) = Ep, (A o f ) for each continuous function /i : f2* — > M, 2 that is 
the two Dynamical Systems (f2*,T,P m ) and (f2, f, P p ) are isomorphic and so are 
the cr-algebras Th = a-{x,coi, . . . ,u>k\ and Tk = o~-{xo, . . . ,Xk}- We will use the 
two processes above interchangeably as far as the study of measure theoretical 
properties is concerned. 

For each function / G C°°(T 2 ,R), / ^ 0, such that rn(/) = we can then define 
the random variables Xk(uj,x) := f(ir o F k (cu,x)), where n(to,x) := x. We are 
interested in studying the P p -almost sure asymptotic behavior, as N — > oo, of the 
random variables 

N-l 

S N (cu) := 5^X fc (w,-). 

fe=0 

The first relevant fact lies in the following Lemma. 
Lemma 2.1. The dynamical system (Cl,F,Pp) is ergodic. 

Proof. By the above discussion the ergodicity of (Cl, F,~P p ) is equivalent to the 
ergodicity of the stationary Markov process P m . It is well known that the ergodicity 
of such a process is equivalent to the fact that Q p g = g implies g = constant for 
each bounded measurable g. In section 3 we will see (Corollary 3.3) that there 
exists p, q > such that, for each / G C p+2d+1 ,g G C q , holds 

(2.1) lim m(fQ n p g) = m(/)m(. 9 ). 

Taking / G C p+2d+1 and g G L°°, we can choose a sequence (gj)j^i C C q that 
converges to g in L 1 and we obtain (2.1) also for such functions. But this means 

2 Indeed, if h(x ,x 1 , . . .) = g(x n ), then E P (h) = M(Qp9) = E P„ „ (h o On the other 
hand if we have already the equality for functions depending on n variables, we can write 
h(xo, ■ ■ . , x n , x n +i) = gx ,...,x„(x n +i) and, by induction, 

EpJh)=E P JE P Jh | xi,...,xn))=K Pll (Q p g ai ,..., Xn (x n )) 

= E Ppi „ x) , . . . , i)„, Ti*(w, x)„) j 

= Ep^ (Bp p ,„ Pk>*I«,«i = Ep p , m (ho*). 

The assertion follows then by the density of the local functions among the continuous ones. 
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that Q p g = g implies m(fg) = m(/)m(g) for each / £ QP+^d+i w \ l [ c \ 1 readily implies 
that g is constant. □ 

Since m(Sjv(w)) = 0, thanks to the previous Lemma, we can apply the Birkhoff 
Ergodic Theorem and obtain 

lim —Sn = P p — almost surely. 

The next step is to investigate the variable N~^Sn and prove that it satisfies a 
(quenched) CLT. 

Theorem 1. For each f £ C°°(T 2 ,R) 7 3 / ^ 0, and for each p £ [0, 1] there exists 
Y? p £ R+ such that, for f p -almost all sequences u> holds 4 

N'^Sn 7V(0, Sp) under m. 

In fact, T, p depends analytically on p. Moreover, if f is not a simultaneous contin- 
uous coboundary for each admissible 6 map Ti, then Y? p > 0. 

Remark 2.2. Imposing Y? p > clearly excludes fewer observables in the random 
case than the deterministic one. From Theorem 1 and the classical Livschitz Theo- 
rem [14] easily follows that Y? p = if and only if^2j— f( x j) = whenever {xj)^~Q 
is a closed orbit for some sequence of admissible maps T Ul , . . . , T Wfc . Unfortunately, 
to use such a criterion it may be necessary to check a very large number of trajecto- 
ries. Yet, if p £ ]0, 1[, then the situation may be much simpler. Indeed, if ' £ 2 = 0, 
then it must be f = g~goTo = g—goT\ and hence goT^oT^ 1 — g. Thus ifTioT^ 1 
is ergodic, 7 then g must be constant and hence f = contrary to assumptions. That 
is: if {To,Ti} are admissible and T\ o T^ 1 is ergodic, then Y? p > 0. 

Remark 2.3. Note that one cannot possibly extend our results to include all se- 
quences; for instance, there exist sequences containing alternating, "deterministic", 
stretches of either Tq 's or T\ 's. If these stretches are of rapidly and ever increasing 
length, then the variance fails to exist. 

Before proving such a strong result we will obtain its averaged (annealed) version. 

Lemma 2.4. For each p £ [0, 1] there exists T, p £ R + such that it holds 

N^^Sn =*> 7V(0, E£) under P p . 

3 In the proof we use / S C r (T 2 , R) for r large enough. Yet, since our bounds for r are far from 
optimal (nor do we strive to optimize them) we see no point in giving an explicit bound for r. 

By A/"(0, E 2 ) we mean the centered Gaussian random variable with variance E 2 . The symbol 
=> stands for convergence in distribution. As usual, Af(0, 0) stands for the measure concentrated 
at zero. 

^By a simultaneous continuos coboundary for a set of maps we mean that there exists a 

continuous function g such that / = g — g o T; , for each map {T{ } . 

6 The map Ti is admissible if it appears with nonzero probability with respect to P p . In our 
case the admissible maps are {Tq,Ti} unless p £ {0, 1}. 

^Note that this may easily fail even if To ^ Ti . Indeed, consider the case 

Ao =Q 3) Al = (l 2) 

Then Ti (Tq _1 (asi, 22)) = (2:1, 22 — xi). In fact, all the functions 3(2:1,0:2) = g(xi) are invariant. 
The identity / = g-goT = g-goT x yields *£%=o fonoF k (uj, x) = g{x)-goT i oT L , N o- ■ •oT Wl (x) 
for i G {0, 1} and hence E| = if g is in L 2 . 
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In turn such a result is based on a fine understanding of the dynamical properties 
of certain transfer operators associated to the process. 

Remark 2.5. Note that one could obtain similar results for any finite collection of 
smooth symplectic hyperbolic maps in any dimension or piecewise smooth maps in 
dimension two. This can be achieved at the price of using in the following section 
the functional setting of [13, 5] or [7] for the piecewise smooth case. 

Our first task will be to obtain some information on the spectral properties of 
such operators. To do so in a useful way it is necessary to introduce appropriate 
functional spaces. Instead of appealing to the general theory developed in [6, 4, 12, 
5, 13] we will take advantage of the simplicity of the present setting and introduce 
explicitly a particularly simple version of such a theory. We will then see how it 
can be used to address the ergodic theoretical questions we are interested in. 

3. Spectral properties of the Transfer operators 

For further use (see section 5) we need to study more general automorphisms 
than the one introduced in the previous section, namely (J)^ T, : T 2d — > T 2d 
defined by (J)^ TiX = (J)^ Aj mod 1 where @f A, is the c?-fold direct sums of the 



matrices Ao,Ai <G SX(2,Z), that is 



e SX(2d,Z), 



1 V Aj 

with d € {1,2}. These matrices are symplectic with respect to the symplectic form 



J Id 



J: 



2, 



where J 2 = ( _?i J) is the standard symplectic form in two dimensions. 8 

Let us introduce the transfer operators C)p (p :— tp o T,) -1 induced by the 
above automorphisms. We will consider the operator obtained from the latter by 
averaging over the Bernoulli measure, namely 

Cf :=p4f + (l"P)4f- 
We also need to study perturbed operators of the form 

£<>:= 4?>(eV)- 
In order to avoid unnecessary proliferation of indices, we set 

(3.1) £r 4 :=4?» C P :=C$>, and C g , p := . 

Finally, notice that the transfer operator C p and the Markov operator Q p are dual: 

(3.2) m{fQ p g) = m(C p f ■ g). 



^In fact all the following can be easily generalized to any set of symplectic toral automorphisms 
which preserves the standard sector in any dimension, see [16, 17] for more details on the necessary 
machinery. 
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To study such operators it is necessary to introduce appropriate Banach spaces 
(see [13, 5, 12]). Here, given the simplicity of the situation, we can quickly introduce 
and use spaces inspired by [5, 12] whereby making the presentation self-consistent. 9 

Given v := (v±, . . . ,V2d) G M. 2d , let us denote v := (ui,t>3, . . . , i>2<z-i) S and 
v '■= i v 2, V4, . . . , V2d)- Then the standard inner product in M. d of the latter two reads 
(f),v) = J2k=o w 2fc+i«2(fc+i)- Consider now the cones C_ := {v G R 2rf : (v,v) < 0}, 
C+ := {v G R 2d : (v,v) > 0}. Then there exists 1 < A < A and Co, depending 
only on {^4o, A\}, such that, for all v G C-, (ii, . . . , i n ) G {0, 1}™, and n G N, 

(3.3) C A < — , — <C A , 

and, for all v G C+, (ii, . . . , i n ) G {0, 1}™, and n G N, 

(3.4) C A < — , — <C A. 

Moreover one can compute that there exists (3 > 1 such that (0^ ^4~ 1 )C_ C Cp, 
where Cp := {w G R 2d : P^M 2 < -(v,v) < f3\\v\\ 2 } C int C_. 

Now, we proceed to define the norm for the Banach space we want to consider. 
Notice that the natural objects in these cones are not vectors but Lagrangian sub- 
spaces. Recall that, given a symplectic form J, a Lagrangian subspace E C M. 2d is 
a d-dimensional subspace such that (v, Jw) = for all v,w G E. 

For our choice of symplectic form, every Lagrangian subspace can also be written 
as the set E = {v : v = — Uv} for a specific symmetric d x d matrix U. Our 
convention is to write a minus sign in front of the U here, because then E C Cp if 
and only if <U </3t. 

Let us denote the set of Lagrangian subspaccs as L. For a Lagrangian subspace 
E and a vector k, we set (E,k) := sup (v,k). Then, for each p,q G M+ and 

IMI=i 

/ G C°°(R 2<i ,R) we define the norm 

H/IU == -P £ IM 1 + |( ^ )r , +l/ol, 
ECC- kez™\{0} lx ' 71 

where fk are the Fourier coefficients of /. 

Notice that the nth power of 4 can be expanded 

n 1 

[4 d) r = E E p 5<i ' o+ " +5i -°(i - p) 5 -^-^- 1 4? . . .4* . 



3=1 i 3 =0 

Because the Bernoulli weights above sum to unity, 



H[4 d) r7iu< sup n4?> ...41! /Hp,, 



(ii,...,i„)e{0,l}" 



9 A ctually our choice is more flexible than the one in [12], in the spirit of [5], and would allow 
to treat C k maps, although it is not the goal here. 
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A straightforward computation shows that f)k = /(® d a t ) k f° r ^ ^ ©1 ^ 2 — 
l? d . Hence, for each n G N and (ii, . . . , i n ) € {0, 1}™, 

I, .(d) .(d) ,n _ ,, , 1(01^ ■■■0i^D~ lfc l P ,,, 

- S i£ 5 {01 l/ ' l i + l ( B ,(e^s...e^ D -., )r « +IA| - 

We begin estimating the summand. Before that, we simplify notation and rename 
the matrix product. To avoid the problem of too many indices, we denote it only 
with the subscript n and assume the dimension and the sequence to be implicit: 

d d 

%i '■= • ■ 1 @ A? 



1 1 
Using (3.3) and the fact E C C_, we have 



(3.5) C A n (E n ,k)>(E,%- 1 k)>Co 1 \ n (E n ,k) where E n := {%~ 1 ) T E 
and consequently, 

I r I \jjKj 1 ^I P < r^p+q I f I \j^n 1 ^I P 

1 'l + 1(^,21,7^)1^+9" Ukl l + X n (P+i)\(E n ,k)\P+i' 

Now there are two possible cases depending on where Sl" 1 ^ lies. If ^ C_, 

then 121-^1 < C AA- n+1 |fc| and we have 

l^fw k\ P pp/ a \ — 1 \p \—np 

l + \n(P+i)\(E n ,k)\P+i ~ ol j 1 + |<£7„,fe)|P+9" 

While if Strife G C_, then 

(3 ' 6) ( ^' fc> -2^yiTF-^, n - 

Indeed, setting k n := 2t~_ 1 fe, 

(v,k n ) - (v,Uk n ) 



((©A rl )^.U > inf = inf sup 



Jet 1 p-ii<u<m vm" WW + WW 

. . „ (k n >k n ) + (k n ,Uk n ) 
> mi — 



where we have chosen v = —k n in the second line. On the other hand the choice 
v = [7 _1 A; rl yields 

((®A£)E,kn)> inf (E,k n ) > 



The inequality (3.6) follows from the above estimates, (3.5) and CoA™ 1 |fc„| > \k\. 
Next, we consider two subcases. If |fc| > Bp^n, then \ (E n , k}\ > 1. Hence, 

j^-ifcjP AP"A-(p+«)"|/c|p AP"A-( p+,? )"|fc|P 
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which is a good estimate provided q is large enough so that A p A _ ^ p+9 - ) < 1. The 
remainder is a finite sum which can be estimated because if E C C_, then E n C 
(%~ 1 ) T C- C C_. Thus, we have 

V lA-M^W < CPA" sup V l/fc||fc|P 

ecc- fcez 2d \M ecc_ fcez 2d \{o} 

Accordingly, setting fi :— max{A _p , A p \~ p ~ q } we can collect all the above inequal- 
ities as 

iii^r/n,,. < c,n/iip„, 
(3 . 7) ii[^>]"/ii„<^vii/ii,„ +s y; T+ifZy +IAI 

|k|<B/9,n 

< C 2 /i n ||/||p, 9 + B n ||/||p_l, 9+ l, 

where B n = C^ +2p Bp_ n A np . Next, for each /i e (jl, 1) choose n such that C^/i" < 
and, for each n € N, write n = kno + m with m £ {0, . . . ,uq — 1}. One can 
thus iterate the second of the (3.7) and obtain 

fc-i 

||[4 d) ]"/IU < M fcno ll>C m /|| P . 9 + B„ EM^II^-'-^^/llp-Lg+i, 

j'=o 

which finally yields 

ii[4'°] n /iu^ c, iii/iu 

||[4 d) ]"/!lp, 9 < C '3M' l ||/ll M + S||/llp-i, (;+ i, 

with B = C2B no (\ — /i™ ) -1 . We can then consider the closure, B p,q , of C°° in the 
space of distributions with respect to the norms |j • \\ p , q . ft is easy to prove the 
following: 

Lemma 3.1. The operators are well defined bounded operators on B p,q , pro- 
vided A p < X p+q . In addition, the unit ball ofB p ' q is relatively compact in B p ~ 1,q+1 . 

Theorem 3.2. // A p < X p+q , the operator C p acting on B p q has an essential 
spectral radius smaller than fj,. The rest of the spectrum consists of finitely many 
eigenvalues of finite multiplicity, all in the unit disk. The only eigenvalue of modulus 
one is one and the constant function equal to one is the corresponding eigen] 'unction. 

Proof. Lemma 3.1, the Lasota-Yorke type inequalities (3.8) imply the result (see 
[3, 6] for more details). □ 

Before proceeding, let us mention that for any r, n £ N we endow the space 
C r (T\M.) with the norm \\g\\ C r := £l=o Wg^lU 

Moreover, a simple computation shows that C r C B p,q , provided r > p + 2d. 



Corollary 3.3. The equation (2.1) holds true. 
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Proof. First of all notice that, for all / G B pq and g G C q holds 

|m(/ fl )|<X;i/illfl-il<IMIc« [I/0I+ e i/hi^i 

V iSZ 2 \{0} 



< 



IpIIc^H/IIp.q max 



fi — W+ q — 1 ^ 2 ^ c "^ q - 

\ ;gz 2 \{o} 



If p, q satisfy the hypothesis of Theorem 3.2, C p has a spectral gap 8 p > 0. Thus, 
|m(/Q» - m(/)m(<?)| = |m(££/ • 5) - m(/)m(fl)| < C(l - 8 P r\\f\\ p Jg\\c; 

because decomposing C p := Q + 1Z with Qf := m(/) we have Q1Z = 1ZQ = and 
therefore £"/ = ft"/ + Qf, where p"^^™) < C(l - S p ) n by the spectral radius 
formula. □ 

Here is the last fact we need to know about the above functional analytic setting. 

Lemma 3.4. For each function g G C 2p+q+2d+1 (T 2d ,M.) the multiplication operator 
M g defined by M g f = gf, is bounded in B pq by C||(7||c2 P +,+2d+i • 

Proof. We define the norm 10 

|| ff || r := sup 1^1(1 + |Af). 

k£Z 2d 

Clearly g E C r {T 2d ,R) implies || 3 || r < |j 5 |j c ,. The S^-norm of the product then 
reads 

(3.9) ||/<7|U=sup E \fi\\9k-i\ „ |ffl, p+g + £ l/'llg-'l- 
e E cc k,lez M ^1^'"/! ;ez2d 

fe#0 

Let us analyze the second term first: 



£ l/j|IS-l|<IMIr l/0|+ E I/' 
JGZ 2d \ ZeZ 2d \{0} 

The desired bound follows, if r > q, from 

1 \l\ p 



J!" 



1 + K-E.J) 



Now, look at the summand in the first term of (3.9), ignoring the I = case that 
can be taken care of separately. 



\k\ p 

M9k - ll l + \(E,k)\ p + q - h 



\fl\ 



l + \(E,l)\ p + q _ 
\k\ p l + \(E,l)\r+* 1 
1 + \(E, k) \ p+q \1\p l + \k-l\ 



10 We are aware that our choices of norms and the subsequent estimates, are not the optimal 
ones. We are simply trying to simplify the arguments as much as possible even at the expense of 
some, not really relevant, optimality. 
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Notice that \(E,l)\ < \(E, k}\ + \ (E, k — l)\. Thus, on the one hand 



k\P l + \{E,k)\*+* 1 <v \k\ p 



, \(E,k)\P+1 \l\P l + \k-l\ r \l\P(l + \k-l\ r ) 

One the other hand 

| fc |P i + \(E,k-l)\^ l v \k\p 

^ 1 + !(£,*;) l + |A:-?| r ~ ^ |Z| P (1 + \k - l\ r -P- q ) ~ ' 

provided r> 2p + q + 2d. The general term of (KE 1 , fc) | + \(E, k - l)\) p+q is bounded 
similarly, using | (E, k)\ n \ (E, k — l) \P+ q - n < (1 + \(E, fc)| n )(l + \(E, k - l)\P+i- n ) for 
each n = 0, . . . , p + q. □ 

4. Averaged CLT 

To establish Lemma 2.4 it suffices to compute 

lim Ep (e'^v ^™=° 1/ '" 



N- 

where f Ut k '■= Xk(u>, •) := / o T Wfc o • • • o T Wl with /^o : = /, and show that this 
limit is the characteristic function e 2 & of the centered normal distribution with 
some variance E^,. Recalling the transfer operators in (3.1), 

(4.1) E P (e i ^w j: "=° f "A=m(C N , 1). 

Hence Theorem 3.2 and Lemma 3.4 show that £ i is a bounded operator 

i\N-?f, P 

on B p q provided / e £2p+g+5 an( ^ depends analytically on A. To continue it is 
necessary to study the leading eigenvalue of such an operator. 

Note that, in general, given any positive operator C on the spaces B p ' q with 
maximal simple eigenvalue one, with a spectral gap and m(Ctp) = m(tp) for each 
smooth ip, for any smooth complex valued function g we can define the family of 
operators C u (p := £(e U9 tp) and, thanks to Lemma 3.4, the standard perturbation 
theory applies. Thus there exists <p v , fi u , with /io = 1, such that 

Cv4>v = Vv4>u, m.(</>„) = l. 

Differentiating this relation with respect to v and integrating one readily obtains 

(4.2) ^mtse^ + e 1 ^) 

and, setting v = 0, <j)' = (1 - C)' 1 [C(gcj) ) - </> m(.gc/>o)] = J2n=o £"[£(#0) ~ 
0o r n(.9 ( / ) o)]- 11 Finally, differentiating again yields 

00 

(4.3) Mo = m(9 2 M + 2 ^ m(gC n [£(g<f> ) - o m(.#o)])- 

Thus, by standard perturbation theory and in view of Theorem 3.2 we can write 
£-v = i^vQ-v + T^u where = Q v , H V Q V = Q v "R-v = 0, the spectral radius of TZ„ is 
smaller than p < 1 for all \i/\ < v$ for some vq > 0, and \p. v — 1— Mo^~ i/ i o 1/2 | — ^Vl 3 
for some fixed constant C > and \v\ < fo- In addition, Q w is a rank one operator 
of the form 4> v £g> m v where m v belongs to the dual of the space, mo = m, and 
K(l)-l| <C\v\. 



The latter is well defined since (1 — C) 1 is applied on a function from which the eigendi- 
rcction of C corresponding to eigenvalue 1 has been projected out, and because of the spectral 
gap. 
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If we apply the above to the operator £ 1 , v = iXN a (hence q = f, and 
0o = 1), then remembering equation (4.1) it follows that 

(4-4) 

E Pp (e^SLo 1 /-,^ = m (^Q, y l+^l) =e Ann ^)(l + 0(|A|^^)) + O(p Ar ) 

= e -^'A 2 +iV ((|A|/iV-i)3) + 0{m -i + p N^ 

because m(Q„l) = m{4> v m v {l)) = m„(l) and |m(7^1)| < \\K^ \\ c{B v. q) < Cp N . 
Hence, if \X\/VN is sufficiently small, we have 

,1 + IAI 3 



(4.5) 



< c- 



N 



Lemma 4.1. The quantity Y? G K defined by 



E„ := lim — E P 

P AT-OO AT " 



'N-l 
lk=0 



II 

Mo- 



is always nonnegative and given by 

oo 

(4-6) m(/ 2 ) + 2 ]>>(/££/) 

n=l 

T7ie map [0, 1] — > R+ : p i— > Sj, is analytic. Moreover, if f is not a C a simultaneous 
coboundary for the admissible automorphisms Ti (see Remark 2.2), then £ 2 > 0. 

Proof. A direct computation yields 



'N-l 

E** 

Lfc=0 



JV-1 



^E Pp (X fe 2 ) + 2 £ e p,(^) 

fc=0 CKjKAKAT-l 

iVm(/ 2 ) + 2 ]T m(/Q*" J '/) 

0<j<A;<Ar-l 



AT 



iV-1 



m(/ 2 ) + 2^m(/£™/) 



iV-l 



2^nm(/£^/) 



Using Corollary 3.3, the last sum converges exponentially fast in n. Hence £ 2 
exists and is nonnegative simply because it is the limit of a nonnegative quantity. 
To address this last issue, suppose Y? — 0. Then 





'N-l 




OO 


JV-1 






] 


< 2N J2 " 


h2^n|m(/£^/)| <C 




_fc=0 




n=JV 


71=1 



uniformly in A". This means that the random variables Zn '■= ^2^ = q Xk are 
uniformly bounded in L 2 . By the Banach-Alaoglu Theorem, they form a weak-* 
relatively compact set. We can then extract a subsequence (Nj)°^ 1 such that, for 
each tp G L 2 (n, P p ), 

lim E Pp (<pZ N .)=E Pp (ipY) 
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for some L 2 random variable Y. If we choose ip to be a function of the x only, it 
follows that 

iV,--l 



lim m(tp V Qlf) = m(ipg) 

— knn * ■ 



7 — >OC 

n=0 



where g = K p (Y) G L 2 (T 2 ,m). On the other hand, for each smooth <p, 

/ / Nj-l Nf-l 

m(<K/ - 9 + Q P 9)) =.limmU/-^ Q p f + £ Q p +1 f 



n=0 n=0 



lim m(<pQ%>f)=0. 



J->oo 

- 2 



That is / = g — Q p g. Next, consider the L random variables G n := gono F n and 
Mn+i = Ek=o( x k + G k+ i - G k ) = G n+1 - Go + ELo x k- For each N e N, we 
use Jensen's inequality to get 



C > E r 

(4.7) 



JV-l 

* k 



k=0 



= JQ,p 



tJ ([M N ~G N +g] 2 ) 



> E Pp (M* ) - 2^E Pp (M 2 N )E Pp ([G N ~g] 2 ). 
In fact, the process (M n ) is a martingale, since 

E Pjs (G fe+1 | x,o;i, . . . ,w fe ) = E Pm (.g(a; fc+ i) | = Q p g(x k ) = (Q p g) otto F fe (.T,w). 
Thus, 

JV-l 

e p, (Mn) = E e p p + G *+i - G *] 2 ) 

fc=0 

= JV {p m([/ + g o T - #) + (1 - p)m([/ + 5 o H - g] 2 )} . 

The inequality (4.7) and the boundedncss of Ep p [[Gn — g} 2 ) imply that pm([f + 
go T - g] 2 ) + (1 - p)m([/ + g o Ti - <?] 2 ) = 0, that is / + 5 o T; - 5 = for each 
admissible Tj. 

The continuity of 5 follows from the usual Livschitz rigidity arguments. 12 
In order to prove analyticity of the variance E 2 with respect to p, first notice 
that there is a positive lower bound on the spectral gap S p appearing in the proof of 
Corollary 3.3 in a complex neighborhood of [0, 1]. Thus, the series in (4.6) converges 



12 In fact, in the present simple case one can provide the following direct proof: clearly g = 
(1 — £t;) 1 C'Tif = y~]uL-[ £7"/ f° r an admissible choice of Ti, convergence taking place in the 
II ' lip. 9 norm. Let v u ' a be the unstable and stable vectors of T;, respectively, and <p G C°°. Then 

00 00 00 

|tn(K, V<p)g)\ < £ |m(^K, V£ Ti /))| = £ W^TTV, C k T . V/»| < £ \F ttU^WvW^- 
fc=i fe=i fc=i 

On the other hand, g(x) = J2k=o f ° T i +9° T ™ +1 > and thc 

mixing of Ti (proven exactly as in 

Corollary 3.3) implies 

oo oo 

|m««',V(p)ff)| < J2 |m(K,V^>/oT^)| < IIV/llooA-^ll^lli!. 

fc = fc = 

Taking the sup over {ip G C°° : HvH^i = l}i it follows that Vg 6 which implies g € W 1 ' 2 . 
Hence, by Morrey's inequality, g G C°. 
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uniformly in p. The partial sums are polynomials of p, hence the limit is an 
analytic function of p. □ 

We finish the section with two simple but important results. 
Lemma 4.2. Denoting f2(x,y) := fix) — f(y), the operator C i satisfies 

iXN 5/ a ,p 



1^ e -A 2 Sj 



* ] 1 



<c- 1 + ' A ' 3 



N 



Proof. The argument follows verbatim the previous discussion. Thus to prove the 
Lemma we only need to compute the second derivative of the leading eigenvalue, 
which we still denote \i v , and to show that /1q = 2E?,. In analogy with (4.3), 



M^m 2 (/ 2 2 ) + 2^m 2 {h[Cf] n f 2 ), 



n=l 

r4 TinoT, rr( 2 )i 



where trt2 is the normalized Lebesgue measure on T . Then [C p '] n f2{x,y) = 
C n p f(x) - C n p f(y), m(l) = 1, and m(/) = yield = 2E* . □ 

Lemma 4.3. There exists Lq > such that, for all L € (0, Lq), the following 
estimate holds, 



(4.8) 




fc=0 



Proof. This is an averaged large deviation estimate and can be obtained exactly 
as the averaged CLT was obtained. Although the idea is standard we give here a 
sketch of the proof. For any random variable Y, for each (3 > 0, 

P P ({Y > L}) < E P Jt {Y > L} e^ Y -V) < e-^Ep p (e^). 

Moreover, P p ({|F| > L}) = P p {{Y > L}) + P p {{Y > —L}). Applying such an 
inequality to the present situation we have 

P P ({Y > L}) < e-P L m(£% N -i f>p l). 

We again apply perturbation theory techniques at the beginning of this section to 
estimate the right-hand side. Using (4.4) with v = ftN^ 1 ,g = f, <f) a = 1, we have 

m(££v-i/,^) = m + <l) = e Nh ^(l + OmN' 1 )) + 0{ P N ). 

If wc define the Legcndre transform Ic{L) = sup|j,| <c Lv — In [i v and we call the 
value in which the sup is attained, then choosing (3 = v*N we have 

P P ({Y > L}) < (1 + C\vA)e- N - Ic[L) + 0{ P N ). 
To compute explicitly Ic(L) we expand 

Lv - h\\i v =Lv — t^Vo + £>(v 3 )- 

Minimizing this quadratic expression leads to a value of = -hr and gives (recalling 
Hq = Sp) the estimate, 

P P ({Y>L})<2e ^ 
provided L < Ce where Ce is small. □ 
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5. Quenched CLT 



Now that we have the CLT in average we would like to establish it for a large 
class of sequences. Let be the variance of the average CLT with respect to 
the Bernoulli process with parameter p. We wish to show that for P p almost all 
sequences to we have the CLT with variance 

To this end we start with an L 2 estimate: assuming that Yjv is a sequence of 
random variables such that Y := limjv E p (Y/v) exists and is real, we can compute 



K(\y n 



Y\ 2 ) = E p (|F^| 2 ) - Y 2 + 2Ydi(Y - E P (Y N )). 



Thus, recalling the notation := / o T u 



E, 



(5.1) 



m(e 



o T Ul and the bound (4.5), 
= E P (|m(e*^ i: ^ 1/w '*)| a ) 
- e-^l + O 



1 + IAI 



N 



The first term on the right-hand side can be conveniently reinterpreted by intro- 
ducing a product system. That is, consider the maps T u>k © T Uk : T 4 — > T 4 , which 

are represented by the block matrices ( A Q k ^4° J S S'i(4, N). Clearly they are 

hyperbolic toral automorphisms (although of a higher dimensional torus) which 
leave Lebesgue measure invariant. The stable and unstable directions are two di- 
mensional. In perfect analogy with the averaged CLT one can define f2{x,y) := 



f(x) — f(y) and study the operator C 



(2) 



(see Section 3). A direct computa- 



tion then shows that, calling rri2 the normalized Lebesgue measure on T 4 , 



E, 



(|m(e 

By Lemma 4.2 and by (5.1), 



E, 



m 2 e 



^Ef = r 1 /2o(r„ t { B r 1J jo...o(T ul fi)r„ 1 ) 



)) 



= m 3 ([£ <2 > , 

VL i\N~1f 2 



E, 



m(e l v / iv 



^k=0 J" 



< C 



1 + IAI 3 



By Chebyshev inequality the above estimate implies 



(5.2) 



{|m( 



e v^v 



'}) 



> £ < Ce 



One would then like to prove almost sure convergence by applying a Borcl-Cantelli 
argument but two problems are in the way: on the one hand the sum over N of 
the above bound diverges, on the other hand one wants the limit to hold almost 
surely for all A, that is one has potentially uncountably many sets to deal with. 
Both problems can be dealt with by applying Borcl-Cantelli to subsequences and 
then showing that controlling the limit of such sequences one controls the limit for 
each N and A. First of all, notice that 



(5.3) 



i-Ly*" 1 f > i^^r N - 1 f , I A — Ai| 

e vTv Zjfc=0 J "' fc — e Vn ^^=0 Ju>,h <- j M 



N-l 
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On the other hand, notice that the estimate (4.8) in Lemma 4.3 also implies 



^ N+M-l 



(5.4) 



VI + MN- 1 - 1 



N+M-l 

E 



i 



N+M-l 



k=0 



> £ 



N+M-l 
^ M-l 

77 E 



> 



2 x /iV + A/[\/l + M^- 1 - 1] 



k=0 



eVN 
2AI 



< Ce 



-CNM- 1 e 2 



Next, consider b g (±, 1) and the sets 13 A k := {2 k + [j'2 t ' A; ] } J < 2 (i_b )fc , A fe := {-k + 
jk- 1 }^- and B k := A k x A*. For each (N, A) € B fe let A k (N, A) = {(Ni,Ai) G 
Nxl : |iV - JVi| < 2 bk + 1, | A - Ai| < k' 1 }. Clearly 

|J A k (N,X) D { (N, A) eNxl : 2 k < N < 2 k+ \ |A| < k) =: J k . 

(N,\)eB k 

We can then write 



r p < sup 

\ [(N,X)£J k 

< E 

(N,A)£B h 



m(eVjv ^'=o -W) _ e ~2 A s, 



> 4e 



sup 

(N lt Xi)eA h (N,X) 



> 4e 



- E | p p ({ m ( e '^ E *= o1 

(AT,A)GB fc I 



> e 



}) 



m 



sup 

.(Ni ,Ai)€Afc(A^,A) 



sup 

(JVi,Ai)GA fc (AT,A) 



> e 



where we have assumed 1 < e in order to deal with the difference e 2 A s fJ _ 

g-2 A i E P . For each bounded function g > holds 

P P ({9 > A}) = E Pp (l {s > A} ) > E p (m(l {g > A} )l {m(9) > 2A} ). 

But m(.g) < |5| 00 m({g > A}) + A, and m{g) > 2A implies m({g > A}) > A\g\^. 
Thus, 



,({m(g) > 2A}) < A^g^P^Ug > A}). 



13 Here fxl stands for the integer closest to x. 
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We can estimate the above expression by 



sup 

(N,X)£J k 



1 ,/w ^1=0 



m(e vTv 



< £ 

(N,X)tB k 

+ 4e- 1 P 



e 2 f 



')-e- 



> 4e 



> e 



}) 



sup 

{N u \-L)eA k {N,X) 



sup 

(N 1 ,X 1 )eA k (N,\) 



e 



> 



Thus, remembering (5.3), 



(I 



sup 

(N,X)e.7 k 



m(e v^v ^<=o e - 



> 4e 



* £ 

(N,\)EB k 

+ 4e- 1 P t , 



< E 

(N,x)eB k L 



y^JV-l f i 2 V 2 



> e 



}) 



sup I Ai| 

(N 1 ,X 1 )eA k (N,X) 



sup 



|Ai-A| 



1=0 



m(e^ E ™ ° ) - e - ^ 2 ^ | > ej) 



ZV— 1 

£/* 



/(- JVi-1 JV-1 

|JV 1 -W|<2"- + 1 VI ViVl Z=0 ViV i=0 



> 



4A- 



JV-l 

E ^ 

i=0 



> 



ke 



Then the estimates (4.8), (5.4) and (5.2) imply, for k > £ p £ _1 , 



(I * 



sup 

(N,X)&J k 

<c E 

(N,X)eB k 



< Ck 2 2^ k e- 



m(e 



,/Tv z -i=0 J". 



2 A E p > 4e 



T + lAp 



e-^^-t + 2°' c e 



E 



-1 -CVe" 



|JVi-A r |<2 i > fc + l 



bk~C2 k< - 1 - b '>e 2 k- 



for which it follows that the sum over k is finite. By Borcl-Cantclli it follows that 

/ i-k= y ,JV - 1 f , . i X 2 V 2 

N,x)eJ k m{e ^ U =° tu ' l )-e-* x p 

finitely many times with probability one. That is, for each e > 0, there exists a 



the above events |sup( 



> 4e| happen only 
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random variable N e : fl — > N U {oo}, P p -almost surely finite, such that 



sup 

|A|<log 2 JV 



< e for N > N F 



Here we used the fact that, for each fixed N, |A| < log 2 N implies (N, A) e J[i og atj • 
Let us call fi e the bad set of sequences, involving N e = oo. It is an increasing set 
with decreasing e, such that P p (|J e>0 f2 e ) = lim^o Pp(fie) = 0; the bad set is 
independent of e. 

This concludes the proof and establishes the almost sure CLT where almost sure 
means that, fixing any Bernoulli measure, the set of the sequences for which we do 
not have CLT has zero measure. Note, however, that the limit (more precisely, the 
variance) is not constant but depends on p. This is natural since the deterministic 
limits p = and p = 1 generically have different variances and as p varies, the 
variance should interpolate smoothly between these two extremal values, which 
indeed is confirmed by Lemma 4.1. 
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